Theory for p-Wave Feshbach Molecules 
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We determine the physical properties of p-wave Feshbach molecules in doubly spin-polarized 40 K 
and find excellent agreement with recent experiments. We show that these molecules have a large 
probability Z to be in the closed channel or bare molecular state responsible for the Feshbach 
resonance. In the superfluid state this allows for observation of Rabi oscillations between the molec- 
ular and atomic components of the Bose-Einstein condensed pairs, which contains a characteristic 
signature of the quantum phase transition that occurs as a function of applied magnetic field. 

PACS numbers: 03.75.-b, 67.40.-w, 39.25.+k 



Introduction. - The continuous excitement in the 
field of ultracold atoms is to a large extent caused by 
the ever increasing experimental control over the cre- 
ation and manipulation of degenerate quantum gases. 
Arguably the most important example of such control is 
the use of Feshbach resonances with which the interaction 
between the atoms can be manipulated by simply vary- 
ing the magnetic field. The so-called s-wave Feshbach 
resonance occurs when atoms collide with zero orbital 
angular momentum (£ — 0) and is experimentally most 
easily accessible. It has been widely used in the study of 
degenerate Fermi gases, especially for the crossover be- 
tween a Bardeen-Cooper-Schrieffer (BCS) superfluid and 
a Bosc-Einstein condensate (BEC) of diatomic molecules 

A novel challenge is to obtain superfluidity also with 
pairs in nonzero orbital angular momentum states by us- 
ing p-wave (1 = 1), or maybe even d-wave (£ = 2) Fes- 
hbach resonances. Up to now, p-wave resonances have 
been seen in Fermi gases of 40 K 0] and 6 Li atoms Q. 
Most recently, Gaebler et al. have been able to create 
and detect p-wave Feshbach molecules in a gas of 40 K 
with lifetimes on the order of milliseconds on both the 
BEC and the BCS side of the resonance [9]. The life- 
time on the BEC side, where the energy of the Feshbach 
molecule lies below the open-channel continuum, is pre- 
sumably limited by collisional losses, whereas the lifetime 
on the BCS side is determined by the decay of the molec- 
ular state into the open-channel continuum. However, 
upon entering the superfluid regime on the BCS side, the 
Bose-Einstein condensed pairs are stabilized by the pres- 
ence of a Fermi sea jlO] . Now that their dominant decay 
mechanism is absent, the lifetime of the pairs on the BCS 
side is expected to become significantly enhanced, giving 
hope that p-wave superfluidity is indeed within reach. 

There are many exciting aspects about p-wave super- 
fluidity. First of all, a p-wave superfluid has a more com- 
plex order parameter than a s-wave superfluid, due to 
the various possible projections of its angular momentum 
(m=0,±l), which can give rise to anisotropic superfluids. 
Second, a p-wave superfluid undergoes various classical 
and quantum phase transitions as a function of temper- 



ature and magnetic field 11, 12, 13|, 14, 15, lj| • As an 
example, the evolution from the BEC to the BCS side 
of the Feshbach resonance at zero temperature contains 
a topological quantum phase transition [TtJ , rather than 
the crossover known from the s-wave case. Third, the p- 
wave resonances are intrinsically narrow. As a result, the 
two-channel nature of the Feshbach resonance becomes 
very prominent. This can be quantified as follows. 

The superfluid phase near a Feshbach resonance arises 
from the Bose-Einstein condensation of pairs or so-called 
dressed molecules [l8j], which are characterized by the 
linear superposition l^drcssod) = y/Z\tp m ) + y/l - ^IV'aa), 
where \ipm) is the wave function of the bare molecular 
state in the closed channel, and IV'aa) describes the atom 
pairs in the open channel. In the case of an s-wave res- 
onance, the magnitude of the probability Z can be esti- 
mated at resonance (unitarity) by Z ~ 2V2A/tt?7 [lj|, 
where A is the universal BCS gap parameter at unitarity 
and rf 1 is a measure for the energy width of the reso- 
nance. For the extremely broad s-wave resonance of 6 Li 
at 834 G, this leads at a Fermi energy of £p = 380 nK to 
Z ~ 4 • 1CT 5 . For the s-wave resonance of 40 K at 202 G, 
we obtain at the same Fermi energy Z ~ 5 • 10~ 3 , which 
is much larger but still a rather small number. However, 
we will show below that for a p-wave resonance the value 
of Z at unitarity can become as large as 0.7. 

This shows that for typical p-wave Feshbach reso- 
nances the two-channel nature of the dressed molecules 
is very important. Whereas p-wave Feshbach molecules 
have recently been carefully studied experimentally, the- 
oretically the physical properties of these molecules have 
not received much attention, even though they are very 
interesting in their own right. In particular, important 
quantities such as the probability Z , the binding energy 
and the lifetime have not been determined yet in terms 
of the experimentally known Feshbach parameters, i.e., 
the background p-wave scattering length Ob g , the width 
of the resonance AB, and the magnetic moment differ- 
ence between the closed and the open channel A/i. In 
this Letter, we show that given these parameters, all 
two-body properties can be computed exactly and in ex- 
cellent agreement with experiment. Moreover, we per- 
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form a two-channel many-body calculation to compute 
the frequency of coherent Rabi or Josephson oscillations 
between the closed and open channel components of the 
Bose-Einstein condensed pairs. 

Two-body physics. — The two-body problem near a p- 
wave Feshbach resonance can be treated along the same 
lines as the corresponding problem for the s-wave reso- 
nance We consider two fermionic particles that are 
in the same hyperfine spin state and that interact with 
each other through a coupling with a closed molecular 
channel and through a p-wave background interaction in 
the open channel, both with angular momentum projec- 
tion to. The scattering due to the background interaction 
generates a ladder sum, which can be conveniently incor- 
porated into the atom- molecule coupling and gives 
rise to a dressed coupling g m (k) = (V>m|V r ma |V'k + Vv / 2 = 

<?mfcm/[l + i(fcabg) 3 ], where |^ k + ^) ^ s the scattering state 
with momentum k in the open channel, Vma. the coupling 
potential between open and closed channels, and a^ g is 
directly related to the effective hard-core radius Oh c of the 
background potential by a^g = ethc/3 1 / 3 . The subscript 
to refers to the projection of the angular momentum 
to = 0, ±1, such that ko = k z and k±i — ^f(k x Hk y )/^/2 7 
while the coupling constant g m characterizes the strength 
of atom-molecule interactions. Note that the dressed or 
renormalized coupling <? m (k) is exact at low energies. 

With this dressed atom-molecule coupling we can cal- 
culate the two-body molecular self-energy analytically, 
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where £k = fi 2 k 2 /2m & is the kinetic energy with m a the 
mass of a single atom. The dimensionless energy £ is 
given by £ = m a a^ g z/?i 2 . To obtain this result we have 
renormalized the bare self-energy TiY^^z) by subtract- 
ing TiYjf n {Q). This is convenient because in the equation 
for the bound state energy of the dressed molecule [l8j ]. 
namely E - 8% = hE^(E), we can subtract ft£„(0) 
on both sides, which also renormalizes the bare detun- 
ing to the renormalized detuning 5 m and leads to 
E — 5 rn = hTt m (E). The two-body resonance now in- 
deed occurs when the bound-state energy reaches the 
atomic continuum, i.e., at E = and S m — 0. Therefore, 
the renormalized detuning is the experimentally relevant 
quantity and given by 5 m = A[i(B — B m ), with B the ap- 
plied magnetic field and B m the location of the Feshbach 
resonance. Note that in general the Feshbach resonances 
with different values of to need not be located at the same 
magnetic field [l9| . This is because only the projection of 
the total angular momentum is conserved in the collision 
of two alkali atoms. Also note that by using the dressed 
atom-molecule coupling, the right-hand side of Eq. fTj is 
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FIG. 1: Two-body physics near the m = p-wave resonance 
of 40 K located at B = 198 G. The upper left panel shows 
the probability Z to be in the closed Feshbach channel as a 
function of the magnetic field B. In the lower left and the 
upper right panel the black line shows the binding energy 
£ as a function of magnetic field. For positive detuning the 
molecules are not stable, resulting in a finite lifetime r, shown 
in the lower right panel. The open circles are data from the 
experiments of Gaebler el al. 



finite, so no arbitrary, undetermined cutoff needs to be 
introduced. For small energy we obtain that ?L£ m (z) ~ 
( 5l 2 „TO a /?i 2 )[-TO a z/97r?i 2 |a bg | - (-TO a z/?L 2 ) 3 / 2 /6^. The 
'universal' coefficient of (— TO a z/?i 2 ) 3 / 2 is in agreement 
with Refs. [Hill. 

The binding energy of p-wave Feshbach molecules in a 
gas of spin-polarized 40 K atoms near the to — Feshbach 
resonance at 198 G has been accurately measured as a 
function of magnetic field by Gaebler et al. [9j, where 
they find an almost perfectly linear behavior. From the 
slope, which is given by ZAfi, we obtain Afi = 0.175/ib 
with jiB the Bohr magneton. However, for small detuning 
they report tiny deviations in this linear behavior. In 
Fig. 1 we have calculated the binding energy with the 
use of Eq. ([1]), where we extract the relevant Feshbach 
parameters ab g and go from the collisional phase shift of 
the considered jj-wave resonance obtained in Ref. [l9| . 
This results in ab g = 98ao and go — 33a^ 2 h 2 /m a , with 
do the Bohr radius. We see that the agreement with 
our two-body calculation is excellent. The nonlinearities 
arise from the slight variation of Z with magnetic field. 

Having obtained the self-energy and the equation for 
the binding energy of the molecules, we can also deter- 
mine the two-body wave function renormalization factor 
Z, given by Z = [1 - dfi£ m (E) / 'dE]' 1 with E the solu- 
tion of the equation for the binding energy. The behavior 
of Z as a function of magnetic field is shown in Fig. 1 . We 
indeed observe that for the considered p-wave resonance 
Z is very large at unitarity, namely Z ~ 0.7, which proves 
the two-channel nature of the p-wave molecules near this 
resonance. Since the difference in magnetic moment be- 
tween a dressed molecule and two atoms is ZA/i, we see 
that this difference does not go to zero at a p-wave res- 



3 



onance. As a result, the binding energy of the dressed 
molecule vanishes only as (B—B m ) and not as (B—B m ) 2 , 
which is the well-known behavior from the s-wave case 
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In the two-body case, there are only stable molecules 
for negative detuning and strictly speaking Z ceases to 
have meaning for positive detuning. In the many-body 
case the Feshbach molecules are also stable at the BCS 
side, due to the presence of a Fermi sea, preventing decay 
into the atomic continuum [Io| . It is then possible to de- 
termine Z at the BCS side using many-body techniques, 
which has been done for the s-wave case both theoreti- 
cally [l(| and experimentally Q. In the two-body case 
the p-wave molecules are not stabilized by a Fermi sea 
and the decay rate r m is readily calculated from the self- 
energy by r m = — 2lm[H m (E + iO)]. The lifetime is thus 
T m = l/r m . The calculated lifetimes are shown in Fig. 1, 
leading to very good agreement with the experimentally 
observed lifetimes without any adjustable parameters. 

Internal Josephson oscillations. — Having obtained the 
analytic solution to the two-body p-w&ve Feshbach prob- 
lem, the next step is to incorporate this two-body physics 
into the many-body physics. To this end we use the fol- 
lowing effective grand-canonical Hamiltonian with chem- 
ical potential /x a for the atoms and /u m for the molecules, 
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where creates an atom with momentum k, b^ m k creates 
a bare p-wave molecule with angular momentum projec- 
tion m and momentum k, and V is the volume of the gas. 
This Hamiltonian is valid when different values of m lead 
to resonances at well-separated magnetic fields B m , as is 
experimentally the case for the m = resonance of 40 K. 
If this is not true, then on the right-hand side of Eq. ((2|) 
we also have to perform a sum over m. Furthermore, in 
equilibrium we have that /i m = 2/i a , but for reasons that 
become clear shortly, we enforce this relation only at the 
end of the calculations. For p-wave resonances the total 
scattering length satisfies a 3 = a 3 jg — m & gl n /&iTh 2 8 m . If 
we parametrize this in analogy with the s-wave case as 
a 3 = aj^ g [l — AB, n /(B — B m )], we see that the atom- 
molecule coupling obeys g m = H(6tt Afj.ABma^/ 'm^) 1 / 2 , 
with AB m the magnetic field width of the Feshbach res- 
onance. Near resonance the background interaction may 
be safely neglected in the Hamiltonian, except for its ef- 
fect on the atom- molecule coupling g m (k). The reason 
for the latter is that, although the effective Hamilto- 
nian above is designed to be accurate for low energies, 
we still need high-energy states to properly describe the 
two-body physics near a Feshbach resonance. As we have 
seen, the effect of these states is to renormalize the atom- 
molecule coupling and therefore it is essential to use the 



dressed coupling <7 m (k) in the effective Hamiltonian [Io| . 

In order to discuss internal Josephson oscillations 
21 . l22| , we turn to the functional- integral formalism 
18| in which the creation and annihilation operators be- 
come complex fields that not only depend on momen- 
tum, but also on time t. A sketch of the derivation goes 
as follows. First, we rewrite the atomic and molecular 
fields as a k (t) -> a k (t)e ie ^\ b mM {t) -> b mM {t)e iB ^ 
to explicitly consider the atomic and molecular phase 
fluctuations & a (t) and & m (t), respectively. Substituting 
these expressions in the action, it changes in the fol- 
lowing way: g m (k) — > </ m (k)e*[ 29 *w -e,n W] and [ij — > 
fij — TiOj(t), where the last substitution is due to the 
time-derivatives in the action and the index j runs over 
two possible subscripts, namely atomic and molecular. 
Next, we separate the action into the part S p h contain- 
ing the phase fluctuations, and the part S that doesn't 
contain these fluctuations. This allows us to do per- 
turbation theory in the phase fluctuations by expand- 
ing the exponential containing S p h. Expanding in the 
fluctuations up to second order, we see that the coef- 
ficients in front of the phase fluctuations get averaged 
over the original action 5*. For example, the term linear 
in 9 m (t) becomes / dt ro (t) Ek( i m,k(*) 6 m,k(t)), where 
the average EkCrnkWrnkW) equals the total num- 
ber of bare molecules and is conveniently expressed as 
— 9il(/i a , /i m )/9/z m with 0(/i a ,/i m ) the thermodynamic 
potential of the gas. By re-exponentiating all terms up 
to quadratic order, we obtain the effective action for the 
phase fluctuations, 
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where the Josephson coupling is given by the expression 
J = 2j2 w g m (k)(b* (t)a k (t)a^{t))/W. 

In deriving the quadratic action for the phase fluctua- 
tions we made one approximation, namely a gradient ex- 
pansion of the fluctuations that takes into account static 
number correlations. To calculate the coefficients in Eq. 
(|3|) we make another approximation, namely mean-field 
theory. The mean-field thermodynamic potential for a 
Fermi gas near a p-wave resonance at zero temperature 
is the direct generalization of the s-wave case, giving 
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with ?iw k = ^(£k -A«a) 2 + 4|.g m (k)| 2 7i mc and n mc the 
bare molecular condensate density. The gap equation 
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FIG. 2: Josephson frequency ujj as a function of the detuning 
S. Energies are expressed in terms of the Fermi energy ef = 
/2m a , with n = 2n mc + n a the total density of 
atoms, which is kept fixed during the calculation at n — 3 x 
10 18 m -3 , so that ey/Ti ~ 27r x 4 kHz. In the left inset, the 
number of bare molecules n mc and open-channel atoms n a are 
given as a function of the detuning. In the right inset we have 
the atomic chemical potential vs. detuning. 



follows from differentiation with respect to n mc , yielding 
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which can be seen as the many-body generalization of the 
two-body bound state equation E — 5 m = h"E m (E), with 
E = fi m . Within mean- field theory, the coefficients in Eq. 
([3]) are readily determined. For example, —dfl/dfi rn — 
n-mcV, where n mc is obtained from the gap equation. 
Also, we have that —dSl/dfi a = n a V = ^ k [l — (ffk — 
/x a )/^Wk]/2, with n a the atomic density. Even the sec- 
ond derivatives are given by simple analytic expressions 
in terms of momentum integrals. Finally, the Josephson 
coupling is given by J = 2n mc J2k |ffm(k)| 2 /?iwk- These 
expressions all depend quantitatively on ab g , showing 
that the background interaction is important. 

From the effective action we can then derive the equa- 
tions of motion by writing down the Euler-Lagrange 
equations. The remaining second-order differential equa- 
tions can be solved analytically, and gives rise to two 
types of modes. One has frequency zero and is recog- 
nized as the Goldstone mode due to the spontaneous 
broken U(l) symmetry associated with the conservation 
of the total number of atoms in the gas. The other 
mode has a nonzero frequency ojj and is recognized as 
the collective mode corresponding to internal Joseph- 
son or Rabi oscillations between the molecular and the 
atomic components of the Bose-Einstein condensed pairs. 
In Fig. 2, we show the Josephson oscillation frequency, 
the atomic chemical potential, the atomic density and the 
bare molecular density as a function of the detuning with 
the parameters for the m = Feshbach resonance of 40 K. 
Far below resonance, not visible on the scale of the fig- 
ure, we find that ljj ~ S. However, near resonance, where 



the detuning is no longer the dominant energy scale, we 
observe completely different behavior. In particular, on 
approach of the quantum critical point at /i a = 
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we 

find that the slope of the Josephson frequency develops 
a logarithmic divergency. Since internal Josephson os- 
cillations can be accessed experimentally by performing 
a Ramsey- type experiment [231 ] . this characteristic fea- 
ture in the Josephson frequency is an observable signa- 
ture for the occurrence of the quantum phase transition. 
Although the Josephson oscillations are damped by pair 
breaking, we expect this effect to be suppressed since the 
single-atom continuum is Pauli-blocked by the Fermi sea 
and the Bose-Einstein condensed pairs are protected by 
the (mostly) gapped atomic spectrum. Also note that 
the collective mode of Josephson oscillations is intrinsi- 
cally connected to the two-channel nature of p-wave reso- 
nances and is absent in a single-channel treatment of the 
interactions. Finally, we hope that our calculation will 
inspire further research towards p-wave superfluidity. 
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